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Pressure gradients, shell crossing singularities and acoustic 
oscillations - application to inhomogeneous cosmological 
models 
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ABSTRACT 

Inhomogeneous cosmological models have recently become a very interesting alterna- 
tive to standard cosmology. This is because these models are able to fit cosmological 
observations without the need for dark energy. However, due to inhomogeneity and 
pressure- less matter content, these models can suffer from shell crossing singularities. 
These singularities occur when two shell of dust collide with each other leading to 
infinite values of the density. In this Letter we show that if inhomogeneous pressure 
is included then these singularities can be prevented from occurring over the period 
of structure formation. Thus, a simple incorporation of a gradient of pressure allows 
for more comprehensive studies of inhomogeneous cosmological models and their ap- 
plication to cosmology. 
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1 INTRODUCTION 

There has been a significant amount of recent interest 
in inhomogeneous cosmological models as an alterna- 
tive to the ACDM concordance model [for a review 
see ICelerierl lj2007t )]. These inhomogeneous models are 
able to fit many cosmolo gical observations w i thout 
the need for dark energy IIDabrowski fc Hendrvl 1998; 
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lying these schemes is to use inhomogeneous solutions 
of the Einstein field equations rather than the standard 
Priedmann-Robertson- Walker solutions used in ACDM 
models. The implementation of these models suggests that 
we need to live close to the centre of the Gpc-scale void. 
There have already been several methods suggested as 
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a test for these types of models. They are based on the 
time drift of cosmological redshift l|Uzan. Clarkson fc Ellisl 
spectral disto rt ions of the CMB power spectrum 



l|Caldwell fc Stebbinsl 20081). the kinematic S unyaey - 



Zel'dovich effect ijGarcia-Bellido fc Haugbolld l2008bh . 
future measurements of supernova in t he redshift range of 
0.1-0.4 ijClifton. Ferreira fc Landll2008l ) and future Baryon 
Acou stic Oscillations measurements ijBoleiko fc Wvithel 
2008). Based on current observations, the alternative of a 
Gpc-scale underdensity is indistinguishable from the dark 
energy scenario. 

However, there remain a number of concerns regard- 
ing these inhomogeneous cosmologies. In particular, it is 
known that p ressure-free Lemaitre-Tolman ijLemaitrel Il933l : 
lTolman|[T933 ) models can evolve to form shell crossing sin- 
gularities. This is an additional singularity to the Big Bang 
that occurs when two shells of matter collide with each other, 
leading to infinite values of the density. These singularities 
were first discussed in the context of astrophysical appli- 
cations of gravitational collapse, where it has been shown 
that they can form without being hi dden inside an hor i- 
zon, and are therefore globally naked l|Yodzis et al.l fl973). 
However, they are not considered as real singularities for 
a number of reasons. First of all because they are weak 
in the sense that the spacetime can be extended through 
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the singular i ty llNewmanlll986l: Frauendiener fc Kleinlfl995l : 
lNolarJl200A IPlebanski fc Kras inski 200d ). Secondly, as was 
shown bv Ijoshil ljl993l ). an object sent through the singu- 
larity would not be crushed (i.e. they would not be focused 
onto a surface or a line). Finally, in spherically symmet- 
ric, asymptotically flat Einstein- Vlasov systems, these types 
of sin gularities have been proved not to occur ijRein et al.l 
1995). However, this has not been proved in cosmological 
models which are not asymptotically flat. 

In gravitational collapse scenario's, shell crossings gen- 
erally do not require consideration because the initial condi- 
tions are generally unre alistic. However, this is not the case 
in cosmological models. iBoleiko et al.l l|2005h showed that a 
large class of realistic models of voids exhibit shell cross- 
ing singularities when the galaxy wall surrounding the void 
began to form. In many applications, initial conditions are 
chosen such tha t the subsequent evolu tion does not exhibit 
shell crossings l|Hellabv fc Lakdll985l ). however this is ex- 
tremely inconvenient and prevents the study of a wide class 
of models. 

In this Letter we show that the simple incorporation 
of pressure gradients leads to an elimination of shell cross- 
ings over the time scales involved with structure formation. 
Anomalies in the evolution associated with the formation of 
acoustic oscillations (see Sec.[5|) prevent us from determining 
conclusively whether this simple incorporation of inhomoge- 
neous pressure can permanently prevent the shell crossing 
singularities. The structure of the Letter is as follows; In 
Sec. [2] the Lemaitre model is presented. Sec. [3] gives an ex- 
plicit specification of the models being considered. Sec. [4] 
presents the evolution of both models and shows that suf- 
ficiently large pressure gradients can prevent shell crossing 
singularities. Sec. [5] then discuss the occurrence of acoustic 
oscillations due to non-zero pressure gradient. 



2 THE SPHERICALLY SYMMETRIC 
INHOMOGENEOUS SPACE TIME 
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where e is the energy density, p is the pressure, A is the 
cosmological constant, n — 8tyG/c 4 , ' stands for dt and ' 

stands for d r . 

These equations can be reduced to l|Lemaitrel fl933) 
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where M is defined as 

2M(r, t) = R{r, t) + R(r, t)<3~ A{t ' r) R 2 (r, t) - 
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Although in the lit erature the abov e syste m of equations is 
usually credited to M isner fc Sharp! [l964\ it was Lemaitre 
who first studied them. That is why we will refer to the 
spherically symmetric system with inhomogeneous pressure 
as the Lemaitre model. In the zero pressure case t he system 
will be referred to a s the Lemaitre-Tolman model ijLemaitrd 
ll933l : lTolman1ll934h (as is widely accepted). 

The equations of motion, T a/3 -p = 0, reduce to 
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Integrating eq. (J3} the B function can be expressed as 
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where E(r) is an arbitrary function. For the zero-pressure 
Lemaitre-Tolman model, A' = 0, implying the time coor- 
dinate can be rescaled such that e A = 1 without loss of 
generality. Moreover, in this case e B =R ,2 /(1 + 2E). 



A spherically symmetric metric in co-moving 
and synchronous coordin ates can be written as 
ijPlebahski fc Krasihskill2006l ) 
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where dQ. 2 = dO 2 + sin 6d(j> . The Einstein field equations 
for the spherically symmetric perfect fluid distribution (in 
coordinate components) are: 
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3 MODEL SET-UP 

To check how the pressure influences the evolution of the 
spacetime, we concentrate on two models: the pressure- 
free Lemaitre-Tolman model and the Lemaitre model. Both 
models are specified by the same initial data given at the 
last scattering instant. The background model is chosen to 
be a flat Friedmann model with Slmat = 0.3, Ha = 0.7 and 
Ho — 70 km s -1 Mpc -1 . The radial coordinate is defined as 
areal distance R at the last scattering instant: 



R(r,tis 



(10) 



However, for clarity in further use, the ~ sign is omitted and 
the new radial coordinate is referred to as r. In this Letter 
r is always expressed in kpc. 

In this Letter we are interested in how the gradient of 
pressure influences the evolution of cosmic structure, and 
if it can prevent the formation of shell crossing singulari- 
ties. The most common place where these singularities oc- 
cur in pressure-free inhomogeneous models is at the edge of 
voids. This is because the expansion of space inside under- 
dense regions is larger than inside overdense regions and this 
forces matter to flow from the voids towards the surrounding 
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Figure 1. Density distribution for four different instants obtained within the Lemaitre-Tolman model (dotted line) and Lemaitre model 
(solid line). The time is indicated at the top of each panel. R corresponds to the value of the areal distance at the given time instant. Both 
models begin their evolution at the last scattering instant, and they are both defined by the same initial data. At t = 10 9 y differences 
are still small. Between t = 9.95 x 10 9 y and t = 10 x 10 9 y the shell crossing singularity occurs in the Lemaitre model. At this instant, 
the Lemaitre-Tolman model breaks down, however we still present the profile of a function which before the shell crossing was related 
to the density. After the shell crossing, R! < within some region and this leads to unphysical, negative values of p, which is depicted 
in the third panel. 



galaxy walls, implying concentric shells of matter may col- 
lide. We therefore consider a model of a void surrounded by 
an overdense region. Within voids, due to the lower amount 
of matter than in the homogeneous background, the curva- 
ture of the space is negative, thus the explicit forms of mass 
(M) and curvature (expressed by the function E) are 



M = Mo + < 
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for t < 0.5a, 
for 0.5a < £ sC 1.5a 
for 1.5a < £ sC 2a, 
for £ > 2a, 



where Mo is the mass in the corresponding volume of the 
homogeneous universe [i.e. Mo = (AtvG /3c 2 )pb,i s r s and pb,is 
is the background density at the last scattering instant], 



Mi = 8M 2 a~ 3 e 



-3/2 



Mi = -0.3 kpc, a = 12 kpc. 




for r ^ 0.5b, 
for 0.5b sC r sC 1.5b 
for 1.5b SC r SC 2b, 
for r > 2b, 



where E x = 4E 2 & _2 e _1 , E 2 = -1.1 x 10" 5 , b = 10.9 kpc. It 
should be noted that other models of voids are also possible - 
even ones which do not evolve from initial ra r efacti ons but 
from condensation, cf. iMustapha fc Hellabvl l|200ll ). How- 
ever, this particular void model was chosen because it de- 
velops, as we will show, a shell crossing singularity. As can 
bee seen for r > 24 kp<o the mass distribution as well as the 
curvature is the same as in the homogeneous FLRW models. 
These functions were used as an initial condition specified at 
the last scattering instant. The initial density distribution 
for these models is very close to the form given in the first 
panel of Fig[T] One can see here that the void region extends 



from R ~ 1.5 Mp<fl and is surrounded by the galaxy wall 
which has a density up to twice the value of the void. 

We start the evolution of both the pressure-free 
Lemaitre-Tolman and the Lemaitre models from the same 
profile of mass and curvature distributions. The only dis- 
crepancy between these models is with the equation of state, 
which was chosen to be of a polytropic form 

p = Kp 1+1/n . (11) 

The polytropic index is chosen to be n = 3/2 which is 
the case of a mono-atomic gas. This equation of state is 
a good approximation to describe degenerate star cores, gi- 
ant gaseous planets, or even for rocky planets. Thus, al- 
though realistic conditions within high-density regions in- 
side walls might lead to a more complicated dependence 
of pressure, this simple polytropic equation of state can be 
treated as a good first approximation to the problem consid- 
ered in this Letter. The constant K for the Lemaitre-Tolman 
model, which is pressure-free, is K = and for Lemaitre 
model is chosen to be K = 1.98 x 10 14 m 2 /s 2 . (Sec. g]) and 
K = 1.08 x 10 14 m 2 /s 2 (Sec. 0. These are very high val- 
ues. For comparison the ratio of standard pressure of air 
{Pair = 101.325 kPa) to its density [pair = 1.292 kg/m 3 ) 
at T=0°C is approximately equal to 7.84 x 10 4 m 2 /s 2 . Such 
values were chosen in order to better depict the influence of 
pressure gradients on the evolution of matter. However, even 
if such very stiff equations of state are employed, their im- 
pact on the evolution is visible only when density gradients 
become large. Thus, the incorporation of this gradient of 
pressure mostly affects only regions where the shell crossing 
singularities would occur. 

The algorithm which is used to calculate the evolution 
in the Lemaitre- Tolman model is the s ame as the one used 
and described in iBolejko et al . (2005) and the algorithm 



1 For the current instant this corresponds to areal distance of 
R « 26 Mpc. 
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which is used to calculated the evolution in the Lemaitre 
model is the same as the one used and described in iBoleikol 
l|2006l ). 



4 EVOLUTION 

The evolution of matter in both the Lemaitre and Lemaitre- 
Tolman models is presented in Fig. [T] Both models start 
from the same initial conditions at the last scattering in- 
stant. As can be seen at t = 10 9 y (Fig. [TJi) the differ- 
ence between density distributions in these models is still 
small - this is because the pressure gradients are still small0 
However, with subsequent evolution the density contrast in- 
creases, and hence the pressure gradient also increases. At 
t = 9.95 x 10 9 y (Fig. [TJa) the difference between the evo- 
lution of matter within the models with and without pres- 
sure is significantly large. The ratio of p/pb at the maxi- 
mum in the Lemaitre-Tolman model is almost 40, whereas 
in the Lemaitre model it is around 3.7. At this instant in 
the Lemaitre-Tolman model, for r close to the maximum, 
R' is almost zero. Soon, the shell crossing occurs (R' = 0, 
M' / 0). At t = 10 X 10 9 y (Fig. [TJ;) the density in the 
Lemaitre model is everywhere finite - at the maximum 
p/pb ~ 4.1. However, the Lemaitre-Tolman model is not 
applicable, since it breaks down at the shell crossing. This 
is because after the shell crossing, R' < (see curve LT in 
Fig EJl, implying from equation (J6j) that the density becomes 
negative. At the current instant, the density in the Lemaitre 
model is still finite and free from shell crossings. This shows 
that the sufficiently large pressure gradient can prevent shell 
crossing singularities over the period of structure formation 
and their evolution. However as it will be shown in the next 
section it is not certain if simple inhomogeneous pressure can 
permanently prevent the occurrence of these singularities. 



5 ACOUSTIC OSCILLATIONS 

In this section, the constant K is almost half of the value 
used in the previous section. For smaller values of pressure, 
the density increases faster than previously which leads to 
a larger density contrast, which subsequently leads to larger 
pressure gradients. This causes oscillations which are de- 
picted in Fig. [3] The nature of these oscillations is visible in 
the following equation [this is a rearranged eq. ©J 



~ A R 2 



1 . -2 



(l + 2£)exp (-2 J dt^-_ 



R 



+ pR' 



(12) 



As seen from the above expression, when in a region where 
the gradient of pressure is too large, the expansion rate, 
e~ A R 2 , slows down. When one shell slows, then a shell with 
larger r expands relatively faster and the pressure gradient 



10x10 M y 
L2 at t = 11.74x 10 9 y 
LT at t = 10x10% 




Figure 2. The dependence of R'{r). The curve denoted as LI cor- 
responds to the value of R' obtained at t = 10 X 10 9 y within the 
Lemaitre model presented in Sec. [4] (see Fig. [lb). Curve denoted 
as L2 corresponds to the value of R' obtained at t = 11.74 x 10 9 
y within the Lemaitre model considered in Sec. [5] (see Fig. [3j. 
Curve denoted as LT corresponds to the value of R' obtained at 
t = 10 x 10 9 y within the Lemaitre-Tolman model (see Fig. [lb). 



drops, which leads to change of sign of p (initially p was an 
increasing function from to r max - where f'max correspond 
to the maximum). When p' changes sign, the expansion rate 
increases again. This is the origin of the oscillations. How- 
ever, at one stage these oscillations become so large that 
R can change sign and some shells start to collapse. Within 
the model considered in this section, the first shells of matter 
start to collapse at t — 11.76 x 10 9 y (note that the oscilla- 
tions presented in Fig.[3]are shown up to t = 11.74x 10 9 y) At 
this stage the oscillations are so rapid and of large amplitude 
that numerical modelling of this phenomena becomes very 
difficult. This prevents us from further studying the evolu- 
tion and to check if the gradient of pressure is sufficient to 
prevent the shell crossing. The value of R' &tt = 11.74xl0 9 y 
is presented in Fig. [2] (curve L2). 

We also proceeded with smaller values of the spatial 
and time steps. However, after 7? < the analysis becomes 
too difficult to handle properly within the comoving and 
synchronous gauge ([!]), due to the steep gradients in the den- 
sity. This suggests that a more comprehensive and detailed 
study of shell crossings singularities and acoustic oscilla- 
tions should be performed within other coordinate systems. 
For example the Gautreau coordinates for perfect fluid 
spaceti mes discussed in detail and studied bv lLaskv fc Lunl 
l|2007al lbh can provide a more suitable insight. The 
Gautreau coordinates have already proved their usefullness 
in studying shell crossing si n gulari ti es for dust spacetimes 



Newman jl986Tl; Frauendiener fc Kleinl 
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3 It should be noted that it is the gradient of pressure that mat- 
ters. In the case of pressure being only a function of time, the 
function e B has the same form as in the Lemaitre-Tolman model, 
and the evolution of the system with pressure is very s imilar to the 
evolution without pressure - see Bol ejko et all l|2005l ) for details. 



6 CONCLUSIONS 

In this Letter we studied the shell crossing singularities 
which occur in pressure-free spherically symmetric cosmo- 
logical models. These models have recently become very 
popular because they can fit cosmological observations with- 
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Figure 3. Density distribution for five different instants for a 
model with equation of state p = 1 .08 x 10 14 p 5 / 3 [Pa] (see Sec. 
[5] for details). 

out the need for dark energy. However, because of shell cross- 
ing singularities, which occur when two shells of matter col- 
lide, the full parameter space cannot be considered. We have 
investigated the incorporation of pressure into these models. 
Adding a gradient of pressure implies the equations can no 
longer be solved analytically, although conceptually they are 
no more difficult - see eq. (HJ - and they are still simple to 
compute numerically. We showed that the incorporation of 
a pressure gradient can prevent shell crossing singularities 
occurring over the period of structure formation. This will 
enable the unlimited investigation of these inhomogeneous 
cosmological models without the occurrence of these anoma- 
lous singularities. 

We also showed in this Letter that, in some cases, when 
the density contrast is sufficiently high, the existence of pres- 
sure gradients leads to acoustic oscillations. At this stage the 
numerical investigation becomes troublesome. However, by 
choosing a harder form of the equation of state we can post- 
pone the occurrence of shell crossing singularities so they 
do not occur over the time of structure evolution. It should 
be noted that such harder equations of state require that 
the constant K in ifTTjl be unnaturally large - by several or- 
ders of magnitude larger that the realistic value. This means 
that to fully describe the process of acoustic oscillations we 
need to perform a more comprehensive analysis. In partic- 
ular, an application of more realistic fluids which allows for 
an anisotropic pressure and heat flow is essential. 
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